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Abstract 

In weakly inhomogeneous media, Maxwell equations assume a Dirac-like form that is particularly 
apt for the study of paraxial propagation. Using this form, and via the Foldy-Wouthuysen trans- 
formation technique of the Dirac equation, we study the spin transport of paraxial light beams in 
weakly inhomogeneous media. We derive the Berry effect terms and establish the spin Hall effect 
and the Rytov rotation law for polarized paraxial beam transport. 
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I. INTRODUCTION 



The first consistent analysis of the geometrical optics approximation for electromagnetic 
waves in weakly inhomogeneous isotropic media was reported in a work by Rytov [1]. He 
studied the geometric optics to first approximation in the wave number and gave his famous 
law of the rotation of the electromagnetic wave polarization plane. Geometrical attributes 
of this law were described further in detail by Vladimirskii (2j. The Rytov- Vladimirskii 
rotation of the polarization plane, which was observed experimentally in the 1980's [3|, has 
been shown (see and the references therein) to be an example of the geometric Berry 
phase 5). Another fundamental phenomenon associated with the propagation of light in 
weakly inhomogeneous media was discovered by Zel'dovich and co-workers in the 1990's 
6|, |2j. It was termed the 'optical Magnus effect', according to which, waves of right and left 
circular polarization propagate along different trajectories. 

A decade later, Bliokh and co-workers js iO] developed a modified geometrical optics 
approximation for weakly inhomogeneous isotropic media, which contained the Rytov law 
and the optical Magnus effect as manifestations of Berry effects. In this approach, a Berry 
connection (gauge potential) in momentum space stems from the representation in which 
the wave equation is diagonal. The Berry curvature (gauge field strength) associated with 
this connection corresponds to a magnetic-monopole-like field. The connection results in 
a Berry phase that describes the Rytov rotation of the polarization plane. The curvature 
causes the optical Magnus effect, the deviations of waves of right/left polarization being in 
opposite directions and orthogonal to the direction of propagation. This produces a spin 



(polarization) current directed across the principal motion o 



Berry effects have been derived via other approaches too 
to the anomalous and the spin Hall effect of electrons 13 



11 



the beam. Such geometric 



12J. The situation is similar 



17| and is, therefore, also referred 



to as the optical Hall effect or the spin Hall effect of photons. 



Because of their capacity to carry orbital angu 



light beams have found considerable interest 18 



ar momentum in addition to spin, paraxial 



23] . The spin transport of such beams has 



been studied via the standard paraxial wave equation with emphasis on the orbital- angular- 
momentum Hall effect 



251 ] . Now, in weakly inhomogeneous media, where permittivity 



varies adiabatically with position (the permeability being constant), Maxwell equations as- 
sume a Dirac-like form that is particularly apt for the study of paraxial propagation (see 
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next section). In the present work, we take advantage of this Dirac-like form to study 
paraxial spin transport in weakly inhomogeneous media. Our approach is similar to that 
used for relativistic electrons 26] via the Foldy-Wouthuysen transformation technique which 



diagonalizes the Dirac Hamiltonian. We, thus, derive the Berry effect terms and establish 
the spin Hall effect and the Rytov rotation law for polarized paraxial beam transport. 



II. MAXWELL EQUATIONS IN DIRAC-LIKE FORM 



In a free weakly inhomogeneous medium, where the (relative) permittivity e(x) varies 
adiabatically with position (Ve — > 0), Maxwell equations can be written as [27] 



n 



- a* = -m • w 



(i) 



n(x) = \/e(x) being the refractive index. Here, 

/ 



*(x,t) 



\ 



F z 
F z 

F X +iFy J 

where F(x, t) = n(x)E(x, t) +icB(x, t) is the Reimann-Silberstein vector and 

/ 



V 



(2) 




, My 




M, 




(3) 



/ being the 2x2 unit matrix. The matrices Mj are Hermitian and satisfy the algebra of the 
Dirac matrices, including 

Mf = 1, MiMj = -MjMi = iM k k cyclic). (4) 

For monochromatic waves ^(x, t) = e~ luJt ^o(x), Maxwell equations (pQ) reduce to 

(M • p)^ = nV (5) 

where p = — i/c _1 V is the (dimensionless) 'momentum' operator, k = u/c being the wave 
number, p represents the momentum (unit wave vector) of the wave and satisfies the 
standard commutation relations = ik~ x biy (Without loss of clarity, we shall use the 

same symbol for an operator and its eigenvalue to simplify the notation.) In the adiabatic 



approximation, let us write n(x) = no + £(x), where £(x) is a small perturbation around the 
uniform background value n Q . Introducing the notation Mj_ = (M x , M y , 0), p ± = (p x ,p y , 0), 
etc., equation (jSJ) may be written as the Dirac-like equation 

ik-'d^o = H^ (6) 

where 

H = -n p 0M±_ ■ p ± = -no/3 -(f3 + iM ± ■ (z x P± ) 

having used (jl]) in the last expression. Note, however, that the 'Hamiltonian' H is not Her- 
mitian because of the last term, which we shall denote by O. The above form is particularly 
apt for the paraxial approximation, with z as the paraxial direction and the subscript _L 
signifying the transverse plane. In the paraxial approximation we have k± = ^Jk% + k 
(i.e. p± = |p_|_| 1) and \&o(x) = e lknoz tp('x), where ip depends weakly on z. Thus, neglect- 
ing the paraxial derivative of we note that p z = no, as seen immediately by the application 
of (— ik~ d z ) on ^>o- Equation ([6]) then reduces to the Dirac-like eigenvalue equation 

m = -pz*p. (7) 

The first and the fourth components of ip, to be denoted by ip± respectively (see (J2J)), 
represent the two circularly polarized states. Equation ((7|) is the Dirac-like paraxial wave 
equation for monochromatic paraxial propagation in weakly inhomogeneous media. It will 
be used in our study of the spin transport problem in the following sections. 



III. THE FOLDY-WOUTHUYSEN TRANSFORMATION 

In homogeneous media (£ = 0, H = H ), the above Hamiltonian resembles the free 
Dirac Hamiltonian. Like the Dirac Hamiltonian, Hq has an odd (off-diagonal) part O which 
involves Mi and thus couples the circularly polarized states with the ^-component of the 
electromagnetic field F z (see (j3])). We can use the Foldy-Wouthuysen technique to construct 
a matrix U that diagonalizes H via the similarity transformation U~ l H Q U. In this scheme 
U is given by e 1306 , or 

U(p ± ) = exp(M ± • p ± 6) = cosh(p ± 9) + M± P± sinh(p ± 9) 
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having used (j4j). (U is Hermitian and not unitary as in the Dirac theory.) The similarity 
transformation diagonalizes Hq, provided 



so that U~ 1 HqU = -Efi and 



tanh(2p ± #) = ?± 
n 



U{J> ± ) = n ° + E + M ^^ ( 8 ) 



2E(n + E) 



where E = \Jtlq — p^_. Under this transformation, the paraxial wave equation (J7|) thus 
becomes 

H'lf)' = -prf (9) 

where ip' = U~ lr ip and 

H' = U- l HU = -Ef3 - ^ 1 ( P± )/3C(x)t/( P± ). 

The last term on the left hand side of the above can be written in the momentum represen- 
tation as 

U- 1 (p A _)/3U(P±)U- 1 (p±)C(ik- 1 V p )U(p ± ) = (3U 2 C(ik- l V v + k~ l A L ) 

= ^[n Q /3 + iM ± • (z x P± )]C(x + k- l A ± ) (10) 

where 

A ± (p ± ) = iET 1 (p. L )VpE7(p_ L ) = tU-\p L )Vp ± U( Pl _). 
In deriving fllUI) . we have made use of the fact that U~ l (3 = j3U as well as the identity 

G-\x)f(d x )G(x) = f(d x + G- x d x G). 

Direct calculation via OH]) yields 

fiz x p ± i(M ± - P± ) P± iM ± 
±[P±) 2E(n + E) 2E 2 (n + E) 2E ' 

In paraxial propagation, we have E ~ — p z , as p± <C (1 <)tiq. Collecting results, we 
therefore have, 

if/ = (1 - ^-M ± ■ pjV 

^pz 

H' = -Pn(x + k~ 1 A±) - —M ± ■ (z x p ± )C(x + AT X A ± ) (11) 

Pz 



where 



Pz x p , i(M± ■ p ± )p ± iMj 
A ± — 1 r 



4p2 4p3 2p 2 
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IV. ADIABATIC EVOLUTION: THE BERRY EFFECT TERMS 



As mentioned before, terms involving Mi couple the circularly polarized states ifj± with 
the ^-component of the electromagnetic field (which is, after all, small in view of paraxial 
propagation). Such terms are to be ignored in the adiabatic evolution, which means that 
we can project H' on the polarization subspace, keeping only the corresponding elements 
(elements 11,14,41,44). Equations (TIT]) , therefore, read: 

/ 

H' = -a z n(-x + A; _1 A ± ), ip' = 



where a z is the Pauli matrix and 

_ a z z x p ± _ a z z x p 

The paraxial wave equation ([9]), thus, becomes 

Hip a = ap z ip a 

where a = ± represents the polarization and 

ft(x,p) = n(x + ATVAL), A ± = ^-. (12) 

4pj 

The Hamiltonian H governs the adiabatic evolution of polarized paraxial waves. It contains 
the 'spin-orbit' interaction via <jA±, as a result of which, the position operator acquires 
an anomalous contribution according to x — > r = x + k~ 1 crA±. The physical (observable) 
beam coordinates are no longer the canonical coordinates x = ifc~ 1 Vp; they are now r = 
ifc^Dp, where Vp = Vp — icrA± is the covariant derivative. The spin-orbit interaction, 
thus, introduces a Berry connection (gauge potential) A± in the momentum space. This is 
in perfect analogy with the electromagnetic interaction, with the role of the position and 
momentum interchanged (and helicity replaced by charge, of course) and is a well-known 



feature of spin transport |8HlO 



26 



28| . The physical coordinates are non- commutative: 



[n,rj\ = ik 2 <jE ijk B k 

where 

is the Berry curvature (gauge field strength) associated with the connection. (In the last ex- 



pression, which corresponds to a magnetic-monopole-like field, we have used p = \Jp 2 z + p\ 

Pz-) 



V. SPIN HALL EFFECT AND THE RYTOV LAW 



In geometrical optics, the ray equations can be obtained 29[ via the Hamilton's equations 

P = -V X H, x = VpH 

where dot denotes derivative with respect to a beam parameter, which, here, corresponds 
to the beam length. The canonical variables x = r — k~ 1 aA± and p = /c _1 k are considered 
classical, of course. Along the beam, r represents the position and p the (unit) wave vector. 
For paraxial transport, using (TT2]) . we thus have 

p = -V r n(r), r = rVBxp (13) 

which, of course, reduce to the standard ray equations in the zeroth approximation k~ l — > 0. 
In this ('classical') limit, the spin-orbit interaction vanishes and the right/left circularly 
polarized rays follow the same trajectory, which is a straight line along the paraxial direction. 
In the first ('semiclassical') approximation, however, these rays split due to the effect of the 
Berry curvature according to (fT3|) . The deflections from the paraxial trajectory are given by 

p x dp 

la 2p 3 z 

where C is the beam trajectory in momentum space. The resulting displacements, being 
in opposite directions, are orthogonal to the paraxial direction and produce a spin current 
across the direction of propagation. This is the spin Hall effect for paraxial beam transport. 
The phase change suffered by the beam in the course of its propagation is given by 

(f) = —cut + J k ■ dx. = —cut + J k ■ dv -\- u J A± ■ dp 

where, in the last integral we have used the fact that A± • p = 0. This integral represents the 
geometric Berry phase, which is of opposite signs for the two polarizations. For a linearly 
polarized beam, it therefore leads to the rotation of the polarization plane through the angle 

/ A± ■ dp = I - tan 2 9 dipzij f 9 2 dip 

where 9{<^i 1), tp are the zenith and azimuthal angles of the beam, respectively. This is the 
Rytov law for paraxial beam transport. 



Sr ± = k- l o [ 
Jc 
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